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Abstract 

Historically cosmic rays have always been at the intersection of astrophysics with 
particle physics. This is still and especially true in current days where experimenters 
routinely observe atmospheric showers from particles whose energies reach macro- 
scopic values up to about 50 Joules. This dwarfs energies achieved in the laboratory 
by about eight orders of magnitude in the detector frame and three orders of magni- 
tude in the center of mass. While the existence of these highest energy cosmic rays 
does not necessarily testify physics not yet discovered, their macroscopic energies 
likely links their origin to the most energetic processes in the Universe. Explanations 
range from conventional shock acceleration to particle physics beyond the Standard 
Model and processes taking place at the earliest moments of our Universe. While 
motivation for some of the more exotic scenarios may have diminished by newest 
data, conventional shock acceleration scenarios remain to be challenged by the ap- 
parent isotropy of cosmic ray arrival directions which may not be easy to reconcile 
with a highly structured and magnetized Universe. Fortunately, many new experi- 
mental activities promise a strong increase of statistics at the highest energies and a 
combination with 7— ray and neutrino astrophysics will put strong constraints on all 
these theoretical models. This short review is far from complete and instead presents 
a selection of aspects regarded by the author as interesting and/or promising for 
the future. 



1 Introduction 

High energy cosmic ray (CR) particles are shielded by Earth's atmosphere and 
reveal their existence on the ground only by indirect effects such as ionization 
and showers of secondary charged particles covering areas up to many km 2 
for the highest energy particles. In fact, in 1912 Victor Hess discovered CRs 
by measuring ionization from a balloon [1], and in 1938 Pierre Auger proved 
the existence of extensive air showers (EAS) caused by primary particles with 
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energies above 10 15 eV by simultaneously observing the arrival of secondary 
particles in Geiger counters many meters apart [2]. 

After almost 90 years of CR research, their origin is still an open question, with 
a degree of uncertainty increasing with energy [3]: Only below 100 MeV kinetic 
energy, where the solar wind shields protons coming from outside the solar 
system, the sun must give rise to the observed proton flux. Above that energy 
the CR spectrum exhibits little structure and is approximated by broken power 
laws oc E~^\ At the energy E ~ 4 x 10 15 eV called the "knee", the flux of 
particles per area, time, solid angle, and energy steepens from a power law 
index 7 ~ 2.7 to one of index ~ 3.0. The bulk of the CRs up to at least that 
energy is believed to originate within the Milky Way Galaxy. Above the so 
called "ankle" at E ~ 5 x 10 18 eV, the spectrum flattens again to a power law 
of index 7 ~ 2.8. This latter feature is often interpreted as a cross over from 
a steeper Galactic component, which above the ankle cannot be confined by 
the Galactic magnetic field, to a harder component of extragalactic origin. 

Until the 1950s the energies achieved with experiments at accelerators were 
lagging behind observed CR energies which explains why many elementary 
particles such as the positron, the muon, and the pion were first discovered in 
CRs [4]. Today, where the center of mass (CM) energies observed in collisions 
with atmospheric nuclei reach up to a PeV, we have again a similar situation. 
In addition, CR interactions in the atmosphere predominantly occur in the 
extreme forward direction which allows to probe non-perturbative effects of 
the strong interaction. This is complementary to collider experiments where 
the detectors can only see interactions with significant transverse momentum 
transfer. 

Over the last few years, several giant air showers have been detected both in 
ground detectors measuring the secondary shower particles directly in water 
tanks or scintillation counters [5,6] and in fluorescence telescopes detecting 
the nitrogen emission induced by the shower [7,8]. This confirms the arrival of 
CRs with energies up to a few hundred EeV (1 EeV = 10 18 eV), correspond- 
ing to about 50 Joules. The existence of such ultra-high energy cosmic rays 
(UHECRs) poses a serious challenge for conventional theories of CR origin 
based on acceleration of charged particles in powerful astrophysical objects. 
The question of the origin of these UHECRs is, therefore, currently a subject 
of much theoretical research as well as experimental efforts. We refer to Ref. [9] 
for recent brief reviews, and Ref. [10,11] for detailed reviews. 

The problems encountered in trying to explain UHECRs in terms of "bottom- 
up" acceleration mechanisms have been well-documented in a number of stud- 
ies; see, e.g., Refs. [12-14]. In summary, apart from energy draining interac- 
tions in the source the maximal UHECR energy is limited by the product 
of the accelerator size and the strength of the magnetic field containing the 
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charged particles to be accelerated, similar to the situation in man-made ac- 
celerators such as at CERN. These criteria reveal that it is hard to accelerate 
protons and heavy nuclei up to the energies observed even in the most powerful 
astrophysical objects such as radio galaxies and active galactic nuclei. 




Fig. 1. UHECR flux as measured by the HiRes-I and HiRes-II detectors [8], and 
the AGASA experiment [6]. Also shown is a fit to the data of a superposition of 
a galactic and an extragalactic source component. This figure is from the second 
reference in Ref. [8]. 

In addition, nucleons above ~ 70 EeV lose energy drastically due to photo- 
pion production on the cosmic microwave background (CMB) — the Greisen- 
Zatsepin-Kuzmin (GZK) effect [15] - which limits the distance to possible 
sources to less than ~ 100 Mpc [13]. Heavy nuclei at these energies are photo- 
disintegrated in the CMB within a few Mpc [17]. If the sources are not too 
strongly clustered in our cosmologically local environment, a cut-off in the 
spectrum above ~ 70 EeV is therefore expected [16]. However, currently there 
seems to be a disagreement between the AGASA ground array [6] which de- 
tected about 8 events above 10 20 eV, as opposed to about 2 expected from a 
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cut-off, and the HiRes fluorescence detector [8] which seems consistent with 
a cut-off [20], see Fig. 1. The resolution of this problem may have to await 
the completion of the Pierre Auger project [21] which will combine these two 
complementary detection techniques. 

Adding to the problem, there are no obvious astronomical counterparts within 
~ 100 Mpc of the Earth [18,13]. At the same time, no significant large-scale 
anisotropy has been observed in UHECR arrival directions above ~ 10 18 eV, 
whereas there are strong hints for small-scale clustering: The AGASA exper- 
iment has observed five doublets and one triplet within 2.5° out of 57 events 
above 40 EeV [6]. This clustering has a chance probability of less than 1% in 
the case of an isotropic distribution. 

There are currently two possible explanations of these experimental findings: 
The first one assumes negligible magnetic deflection. In this case most of the 
sources would have to be at cosmological distances which would explain the 
absence of nearby counterparts and the apparent isotropy would indicate that 
many sources contribute to the observed flux, where a subset of especially 
powerful sources would explain the small-scale clustering [19]. This scenario 
predicts the confirmation of a GZK cutoff. The second scenario is most likely 
more realistic and takes into account the likely existence of large scale inter- 
vening magnetic fields with intensity B ~ 0.1 — 1 fiG, correlating with the 
large scale galaxy distribution. In this case magnetic deflection would be con- 
siderable even at the highest energies and the observed UHECR flux could be 
dominated by relatively few sources within about 100 Mpc. Here, large scale 
isotropy could be explained by considerable angular deflection leading to dif- 
fusion up to almost the highest energies and the small scale clustering could 
be due to magnetic lensing. This scenario will be discussed in section 3. 

More speculative ways to explain the UHECR observations consist of circum- 
venting the distance restriction imposed by the GZK effect by postulating 
either new particles or new interactions beyond the Standard Model of par- 
ticle physics (section 4) or a violation of the Lorentz symmetry (section 5). 
These possibilities do not, however, solve the problem of acceleration to and 
beyond the observed energies which has to be solved separately. 

In contrast, in the "top-down" scenarios, which will be discussed in section 6, 
the problem of energetics is trivially solved. Here, the UHECR particles are the 
decay products of some super-massive "X" particles of mass mx 10 20 eV, 
and have energies all the way up to ~ m X - Thus, no acceleration mechanism 
is needed. The massive X particles could be metastable relics of the early 
Universe with lifetimes of order of or above the current age of the Universe 
or could be released from topological defects (TDs) that were produced in 
the early Universe during symmetry-breaking phase transitions envisaged in 
Grand Unified Theories (GUTs). If the X particles themselves or their sources 
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cluster similar to dark matter, the dominant observable UHECR contribution 
would come from the Galactic Halo and absorption would be negligible. This 
option seems strongly constrained and may soon be ruled out. 

Non-astrophysical solutions of the UHECR problem are of course in general 
quite model dependent. In addition, if the existence of the GZK cutoff will be 
confirmed, their motivation surely diminishes in general. On the other hand, 
even in this case UHECR can be used to test and constrain new physics beyond 
the Standard Model, such as new interactions beyond the reach of terrestrial 
accelerators (see section 4), violation of symmetries (see section 5), as well 
as Grand Unification and early Universe cosmology, such as the rate of TD 
and/or massive particle production in inflation (see section 6), at energies 
often inaccessible to accelerator experiments. 

The physics and astrophysics of UHECRs are intimately linked with the emerg- 
ing field of neutrino astronomy (for reviews see Refs. [22]) as well as with the 
already established field of 7— ray astronomy (for reviews see, e.g., Ref. [23]). 
Indeed, all scenarios of UHECR origin, including the top-down models, are 
severely constrained by neutrino and 7— ray observations and limits. In turn, 
this linkage has important consequences for theoretical predictions of fluxes 
of extragalactic neutrinos above about a TeV whose detection is a major goal 
of next-generation neutrino telescopes: If these neutrinos are produced as sec- 
ondaries of protons accelerated in astrophysical sources and if these protons 
are not absorbed in the sources, but rather contribute to the UHECR flux 
observed, then the energy content in the neutrino flux can not be higher 
than the one in UHECRs, leading to the so called Waxman Bahcall bound 
for transparent sources with soft acceleration spectra [24,26]. If one of these 
assumptions does not apply, such as for acceleration sources with injection 
spectra harder than E~ 2 and/or opaque to nucleons, or in the top-down sce- 
narios where X particle decays produce much fewer nucleons than 7— rays and 
neutrinos, the Waxman Bahcall bound does not apply, but the neutrino flux 
is still constrained by the observed diffuse 7— ray flux in the GeV range. 



2 Propagation of Ultra-High Energy Radiation 

Before discussing specific scenarios for the UHECR origin we give a short ac- 
count of the interaction and energy loss processes relevant for the propagation 
of ultra-high energy cosmic and 7— rays, and neutrinos [27]. These processes 
have been implemented in various numerical codes used for computing spectra 
of all these particles [28-30]. In the following we assume a flat Universe with 
a Hubble constant of H — 70 km sec _1 Mpc _1 and a cosmological constant 
Qa = 0.7, as favored by current observations. 
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The relevant nucleoli interactions are pair production by protons (jry& — > 
pe~e + ), photo-production of single or multiple pions (N'-fb — * iV nir, n > 1), 
and neutron decay (n — > pe~~v e ). The nucleon threshold energy for single pion 
production on a background photon of energy e is given by 

E th = m - (m " + m - /2) ^ 6.8 x 10 16 f^V 1 eV. (1) 
e VeV / 

Typical CMB photon energies are e ~ 10~ 3 eV, leading to the GZK "cutoff" at 
a few tens of EeV where the nucleon interaction length drops to about 6 Mpc. 

7-rays and electrons/positrons initiate electromagnetic (EM) cascades on low 
energy radiation fields such as the CMB. The high energy photons undergo 
electron-positron pair production (PP; 77^ — > e~e + ). At energies below ~ 
10 14 eV photons interact mainly with the universal infrared and optical (IR/ O) 
backgrounds, between ~ 10 14 eV and ~ 100 EeV mainly with the CMB, while 
above ~ 100 EeV the main target is the universal radio background (URB). In 
the Klein-Nishina regime, where the CM energy is large compared to the elec- 
tron mass, one of the outgoing particles usually carries most of the initial en- 
ergy. This "leading" electron (positron) in turn can transfer almost all of its en- 
ergy to a background photon via inverse Compton scattering (ICS; e% — > e'7). 
EM cascades are driven by this cycle of PP and ICS. The energy degradation 
of the "leading" particle in this cycle is slow, whereas the total number of 
particles grows exponentially with time. Apart from these interactions of first 
order in the EM coupling, triplet pair production (TPP; ejb — > ee~e + ), and 
double pair production (DPP, 775 — > e~e + e~e + ) can play a significant role in 
the formation of 7-ray spectrum in the energy range 10 8 eV < E < 10 25 eV. 
Finally, electrons are subject to synchrotron losses in large scale extragalactic 
magnetic fields (EGMF), whereas for protons synchrotron losses are negligible 
outside the sources. 

Similarly to photons, ultra-high energy (UHE) neutrinos give rise to neutrino 
cascades in the primordial neutrino background via exchange of W and Z 
bosons [31]. Besides the secondary neutrinos which drive the neutrino cas- 
cade, the W and Z decay products include charged leptons and quarks which 
in turn feed into the EM and hadronic channels. Neutrino interactions be- 
come especially significant if the relic neutrinos have masses m v in the eV 
range and thus constitute hot dark matter, because the Z boson resonance 
then occurs at an UHE neutrino energy E Tes = 4 x 10 21 (eV/m y ) eV. In fact, 
the decay products of this "Z-burst" have been proposed as a significant source 
of UHECRs [32]. The big drawback of this scenario is the need of enormous 
primary neutrino fluxes that cannot be produced by known astrophysical ac- 
celeration sources [29]. Even more exotic top-down type sources such as X 
particles decaying mostly into neutrinos appear ruled out due to a tendency 
to overproduce the diffuse GeV 7— ray flux observed by EGRET [33,30,34]. 
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The two major uncertainties in the particle transport are the intensity and 
spectrum of the URB for which there exist no direct measurements in the 
relevant MHz regime [35,36], and the average value of the EGMF. Simulations 
have been performed for different assumptions on these uncertainties. A strong 
URB tends to suppress the UHE 7-ray flux by direct absorption whereas 
a strong EGMF blocks EM cascading (which otherwise develops efficiently 
especially in a low URB) by synchrotron cooling of the electrons. For the 
IR/O background we used the most recent data [37] for simulations in the 
present paper. 

In top-down scenarios, the particle injection spectrum is generally dominated 
by the "primary" 7-rays and neutrinos over nucleons. These primary 7-rays 
and neutrinos are produced by the decay of the primary pions resulting from 
the hadronization of quarks that come from the decay of the X particles. In 
contrast, in acceleration scenarios the primaries are accelerated protons or 
nuclei, and 7-rays, electrons, and neutrinos are produced as secondaries from 
decaying pions that are in turn produced by the interactions of nucleons with 
the low energy photon background. 



3 Deflection in Cosmic Magnetic Fields and the Anisotropy Prob- 
lem 

Magnetic fields are omnipresent in the Universe, but their origin still lies in the 
dark [38]. Best known are the magnetic fields in galaxies which have strengths 
of a few micro Gauss, but there are also some indications for fields correlated 
with larger structures such as galaxy clusters [39]. Magnetic fields as strong 
as ~ 1/iG in sheets and filaments of the large scale galaxy distribution, such 
as in our Local Supercluster, are compatible with existing upper limits on 
Faraday rotation [39-41]. It is also possible that fossil cocoons of former radio 
galaxies, so called radio ghosts, contribute to the isotropization of UHECR 
arrival directions [42]. 

Contrary to the case of electrons, for charged hadrons deflection is more im- 
portant than synchrotron loss in the EGMF. To get an impression of typical 
deflection angles one can characterize the EGMF by its r.m.s. strength B 
and a coherence length l c . If we neglect energy loss processes for the mo- 
ment, then the r.m.s. deflection angle over a distance r in such a field is 
9(E,r) ~ (2r/ c /9) 1 / 2 /ri / [43], where the Larmor radius of a particle of charge 
Ze and energy E is tl ~ E/(ZeB). In numbers this reads 

r » - °- 8 ° z Go^v) ~\m&) V2 (lib) GiRo) • (2) 
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for r > l c . This expression makes it immediately obvious why a magnetized 
Local Supercluster with fields of fractions of micro Gauss prevents a direct 
assignment of sources in the arrival directions of observed UHECRs; the de- 
flection expected is many tens of degrees even at the highest energies. This goes 
along with a time delay r(E,r) ~ r6(E,df/A ~ 1.5 x 10 3 Z 2 (E/10 20 eV)" 2 
(r/10 Mpc) 2 (/ c /Mpc)(_B/10 -9 G) 2 yr which can be millions of years. A source 
visible in UHECRs today could therefore be optically invisible since many 
models involving, for example, active galaxies as UHECR accelerators, pre- 
dict variability on shorter time scales. 

On the other hand, an EGMF of this size could explain the observed large 
scale UHECR isotropy by diffusion and the small-scale clustering by magnetic 
lensing, even if most of the sources are relatively nearby. In fact, numerical 
simulations of nucleons in a simplified scenario where the UHECR source den- 
sity in the Local Supercluster is idealized as a Gaussian sheet of a few Mpc 
thickness and about 20 Mpc length containing a magnetic field with a Kol- 
mogorov spectrum and an r.m.s. strength B proportional to the same profile 
have lead to the following result: About 10 sources in the Local Supercluster 
and a maximal field strength of ~ 0.3/iG lead to arrival direction multi-pole 
moments and autocorrelation functions consistent with the AGASA data [44]. 

Most recent results indicate that the situation is similar in more realistic 
models of the Local Supercluster. We have used the (50 Mpc 3 box of a large 
scale structure simulation [45] where the magnetic field was followed passively 
and normalized to the known strengths of a few micro Gauss in the center of 
galaxy clusters. In this simulated box we selected as observer a place typical 
for our neighborhood in the Local Supercluster, and chose randomly a certain 
number of sources with probability proportional to the local baryon density. 
For each such configuration 5000 nucleon trajectories originating from the 
sources (assumed to emit equal spectra and power) and reaching the observer 
were computed numerically, taking into account energy losses and deflection in 
the magnetic field in the box. For each realization these 5000 trajectories were 
used to construct arrival direction probability distributions taking into account 
the solid-angle dependent exposure function for the respective experiment and 
folding over the angular resolution. From these distributions mock data sets 
consisting of N observed events were dialed. For each such mock data set or 
for the real data set we then obtained estimators for the spherical harmonic 
coefficients C(l) and the autocorrelation function N(6). The estimator for C(l) 
is defined as 

where Ui is the total experimental exposure at arrival direction u l , J\f = 
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57 events with 40. < E/EeV < 9000. 
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Fig. 2. The angular power spectrum C{1) as a function of multi-pole i, obtained for 
the AGASA exposure function, for N = 57 events observed above 40 EeV, sampled 
from 12 simulated configurations of 100 sources in the simulation box of Ref. [45] , 
emitting a proton spectrum oc E~ 2A up to 10 21 eV. The magnetic field strength at 
the observer is ~ 0.13/iG. The diamonds indicate the average over 12 realizations, 
and the left and right error bars represent the statistical and total (including cosmic 
variance due to different realizations) error, respectively, see text for explanations. 
The histogram represents the AGASA data. The overall likelihood significance is 
~ 0.15. 

J2iLi is the sum of the weights 1/uJi, and Yi m {u l ) is the real- valued spheri- 
cal harmonics function taken at direction u l . The estimator for N(9) is denned 
as 

1 ,__ 1 if 6L- is in same bin as 8 
^\y) o otherwise 




and S(8) is the solid angle size of the corresponding bin. 

The different mock data sets in the various realizations yield the statistical 
distributions of C(l) and N(8). One defines the average over all mock data 
sets and realizations as well as two errors: The smaller error (shown to the left 
of the average in Figs. 2 and 3) is the statistical error, i.e. the fluctuations due 
to the finite number N of observed events, averaged over all realizations, while 
the larger error (shown to the right of the average in Figs. 2 and 3) is the "total 
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Fig. 3. As Fig. 2, but for the angular correlation function N{9) as a function of 
angular distance 9, using a bin size of A9 = 1°. The overall likelihood significance 
is ~ 0.63. 



error", i.e. the statistical error plus the cosmic variance, in other words, the 
fluctuations due to finite number of events and the variation between different 
realizations of observer and source positions. For a given data set a x 2 summed 
over all bins can be defined relative to the mean and variance of C(l) and N(9) 
simulated for a given model. From this one can obtain an overall likelihood 
for the consistency of a given real data set with a given model by counting the 
fraction of simulated mock data sets with a \ 2 larger than for the real data. 

We find that as long as the observer is surrounded by magnetic fields above 
about 0.1/iG, which is possible but not obvious for our environment, roughly 10 
or more sources lead to multi-poles and autocorrelations marginally consistent 
with present data which are limited to the Northern hemisphere, although con- 
sistency of large scale multi-poles is somewhat worse than for more extended 
EGMFs. An example for 100 sources is shown in Figs. 2 and 3 [46]. A proton 
injection spectrum roughly oc E~ 2A extending up to ~ 10 21 eV can reproduce 
the sub-GZK spectrum and tends to predict a spectrum somewhere between 
the AGASA and HiRes observations above GZK energies, see Fig. 4. In con- 
trast, if the observer is in a region of EGMF field strength much smaller than 
~ 0.1/iG, the UHECR sky distribution reflects the highly structured large 
scale galaxy distribution and is thus inconsistent with the observed isotropy. 
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Fig. 4. Predicted spectrum observable by AGASA for the scenario discussed in the 
text and shown in Figs. 2 and 3, averaged over 12 realizations, as compared to the 
AGASA (dots) and HiRes-I (stars) data. 

We also find that sources outside our Local Supercluster do not contribute sig- 
nificantly to the observable flux if the observer is immersed in magnetic fields 
above about 0.1/xG and if sources reside in magnetized clusters and super- 
clusters: For particles above the GZK cutoff this is because sources outside 
the Local Supercluster are beyond the GZK distance. Sub-GZK particles are 
mainly contained in the magnetized environment of their sources and thus 
exhibit a much higher local over-density than their sources. Further, the sup- 
pressed flux of low energy particles leaving their environment is largely kept 
out of our Local Supercluster by its own magnetic field [44]. Both effects can 
be understood qualitatively by matching the flux j(E) in the unmagnetized 
region with the diffusive flux —D(E)Vn(E,r) in terms of the diffusion con- 
stant D(E) and the density n(E,r) of particles of energy E which shows that 
the density gradient always points to the source. A significant contribution 
from sources at cosmological distance would require magnetic fields in the 
nano Gauss range, and even then would be unlikely to explain the observed 
clustering. 

The confidence levels that can be obtained with this method for specific mod- 
els of our local magnetic and UHECR source neighborhood will greatly in- 
crease with the increase of data anticipated from future experiments. Full sky 
coverage alone will play an important role in this context as many scenarios 
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Fig. 5. The UHECR arrival direction distribution above 40 EeV in terrestrial co- 
ordinates for the scenario of Figs. 2-4, averaged over all 12 realizations of 5000 
trajectories each. The color scale represents the integral flux per solid angle. The 
pixel size is 1° and the image has been convolved to an angular resolution of 2.4° 
corresponding to the approximate AGASA angular resolution. 

predict large dipoles for the UHECR distribution. This is the case for the spe- 
cific simulation box on which Figs. 2-4 are based, as demonstrated in Fig. 5. 
Whereas current northern hemisphere data are consistent with this scenario 
at the ~ 1.5sigma level, a comparable or larger exposure in the southern 
hemisphere would be sufficient in this case to find a dipole at several sigma 
confidence level. 

Modeling our cosmic neighborhood and simulating UHECR propagation in 
this environment will therefore become more and more important in the com- 
ing years. This will also have to include the effects of the Galactic magnetic 
field and an extension to a possible heavy component of nuclei. For first steps 
in this direction see, e.g. Ref. [47] and Ref. [48], respectively. 



4 New Primary Particles and Interactions 

A possible way around the problem of missing counterparts within accelera- 
tion scenarios is to propose primary particles whose range is not limited by 
interactions with the CMB. Within the Standard Model the only candidate 
is the neutrino, whereas in extensions of the Standard Model one could think 
of new neutrals such as axions or stable supersymmetric elementary particles. 
Such options are mostly ruled out by the tension between the necessity of a 
small EM coupling to avoid the GZK cutoff and a large hadronic coupling to 
ensure normal air showers [49] . Also suggested have been new neutral hadronic 
bound states of light gluinos with quarks and gluons, so-called R-hadrons that 
are heavier than nucleons, and therefore have a higher GZK threshold [50] , as 
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can be seen from Eq. (1). Since this too seems to be disfavored by accelerator 
constraints [51] we will here focus on neutrinos. 

In both the neutrino and new neutral stable particle scenario the particle 
propagating over extragalactic distances would have to be produced as a sec- 
ondary in interactions of a primary proton that is accelerated in a powerful 
active galactic nucleus which can, in contrast to the case of EAS induced by 
nucleons, nuclei, or 7— rays, be located at high redshift. Consequently, these 
scenarios predict a correlation between primary arrival directions and high 
redshift sources. In fact, possible evidence for a correlation of UHECR arrival 
directions with compact radio quasars and BL-Lac objects, some of them 
possibly too far away to be consistent with the GZK effect, was recently re- 
ported [52]. The main challenge in these correlation studies is the choice of 
physically meaningful source selection criteria and the avoidance of a poste- 
riori statistical effects. However, a moderate increase in the observed number 
of events will most likely confirm or rule out the correlation hypothesis. Note, 
however, that these scenarios require the primary proton to be accelerated up 
to at least 10 21 eV, demanding a very powerful astrophysical accelerator. 

4-1 New Neutrino Interactions 

Neutrino primaries have the advantage of being well established particles. 
However, within the Standard Model their interaction cross section with nu- 
cleons, whose charged current part can be parametrized by [53] 

a s J{E) ~ 2.36 x 1(T 32 (£/10 19 eV) 0363 cm 2 , (5) 

for 10 16 eV < E < 10 21 eV, falls short by about five orders of magnitude to 
produce air showers as they are observed. However, it has been suggested 
that the neutrino- nucleon cross section, o~ u n, can be enhanced by new physics 
beyond the electroweak scale in the CM frame, or above about a PeV in the 
nucleon rest frame. Neutrino induced air showers may therefore rather directly 
probe new physics beyond the electroweak scale. 

One possibility consists of a large increase in the number of degrees of freedom 
above the electroweak scale [54] . A specific instance of this idea appears in the- 
ories with n additional large compact dimensions and a quantum gravity scale 
M4 +n ~ TeV that has recently received much attention in the literature [55] 
because it provides an alternative solution to the hierarchy problem in grand 
unifications of gauge interactions without a need of supersymmetry. One of 
the largest contributions to the neutrino-nucleon cross section turns out to 
be the production on the brane representing our world of microscopic black 
holes which extend into the extra dimensions. The production of compact 
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branes, completely wrapped around the extra dimensions, may provide even 
larger contributions [56]. The cross sections can be larger than in the Standard 
model one by up to a factor ~ 100 for reasonable parameters [57]. However, 
this is not sufficient to explain the observed UHECR events [58]. 
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Fig. 6. [59] Cosmogenic neutrino flux per flavor (thick line, assuming maximal mix- 
ing among all flavors) produced by primary proton flux from the AGASA [6] and 
HiRes [8] experiments above 3 x 10 18 eV (error bars). The UHECR sources were as- 
sumed to inject a proton spectrum oc E~ 2 up to 10 22 eV with luminosity oc (1+z) 3 up 
to z = 2. For comparison shown are the atmospheric neutrino flux [60], as well as ex- 
isting upper limits on the diffuse neutrino fluxes from MACRO [61], AMANDA [62], 
BAIKAL [63], AGASA [64], the Fly's Eye [65] and RICE [66] experiments, and the 
limit obtained with the Goldstone radio telescope (GLUE) [67], as indicated. The- 
oretical neutrino flux upper limits based on the 7— ray bound derived from the 
EGRET diffuse 7-ray flux (shown on left part between 30 MeV and 100 GeV) [33], 
the MPR limit for optically thin sources [26] , and the WB limit for AGN-like redsfhft 
evolution [24] are also shown. 



However, the UHECR data can be used to put constraints on cross sections 
satisfying cr vN (E > 10 19 eV) < 10 _27 cm 2 . Particles with such cross sections 
would give rise to horizontal air showers which have not yet been observed. 
Resulting upper limits on their fluxes assuming the Standard Model cross sec- 
tion Eq. (5) are shown in Fig. 6. Comparison with the "cosmogenic" neutrino 
flux produced by UHECRs interacting with the CMB then results in upper 
limits on the cross section which are about a factor 1000 larger than Eq. (5) 
in the energy range between ~ 10 17 eV and ~ 10 19 eV [68-70]. The projected 
sensitivity of future experiments shown in Fig. 7 below indicate that these 
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limits could be lowered down to the Standard Model one [70]. In case of a 
detection of penetrating events the degeneracy of the cross section with the 
unknown neutrino flux could be broken by comparing the rates of horizontal 
air showers with the ones of Earth skimming events [71]. 

The CM energy reached by an UHECR nucleon of energy E interacting with an 
atmospheric nucleon at rest is \fs ~ 0.4(i?/10 20 eV) PeV. Since the starting 
points in the atmosphere of observed EAS at these energies are consistent 
with normal hadronic cross sections, scenarios predicting dramatic increases 
or saturation of cross sections in the TeV range can probably be excluded. It 
is generally expected that new physics appears at TeV scales, and dramatic 
effects are possible when this scale is related to a true fundamental gravity 
scale. There are, for example, string-inspired scenarios with infinite-volume 
extra dimensions and a string scale M s ~ M e 2 ff /M P i ~ (M eff /TeU) 2 10~ 3 eV, 
where M e g is the energy scale where cross sections would undergo a Hagedorn- 
like saturation [72,73], and Mpi ~ 1.72 x 10 18 GeV is the reduced Planck 
mass. The UHECR argument would then require M e g > PeV and thus M s > 
100 eV. However, in case of at least two infinite extra dimensions gravity would 
appear higher-dimensional, i.e. decreasing with a higher power of distance than 
Newtonian gravity on scales r > M P1 /M S 2 ~ 1 pc (M s /100eV)~ 2 . The above 
bound would thus imply modifications of gravity on parsec scales which seems 
phenomenologically inviable. This is an example where UHECR observations 
can constrain models not constrained yet by collider physics due to the lower 
CM energies reached. Of course, UHECR observations are much more indirect 
and thus are only sensitive to dramatic effects such as Hagedorn saturation 
and macroscopic cross sections [74]. 



5 Violation of Lorentz Invariance 

The most elegant solution to the problem of apparently missing nearby sources 
of UHECRs and for their putative correlation with high redshift sources would 
be to speculate that the GZK effect does not exist theoretically. A number of 
authors pointed out [75,76] that this may be possible by allowing violation of 
Lorentz invariance (VLI) by a tiny amount that is consistent with all current 
experiments. At a purely theoretical level, several quantum gravity models 
including some based on string theories do in fact predict non-trivial modifi- 
cations of space-time symmetries that also imply VLI at extremely short dis- 
tances (or equivalently at extremely high energies); see e.g., Ref. [77] and ref- 
erences therein. These theories are, however, not yet in forms definite enough 
to allow precise quantitative predictions of the exact form of the possible VLI. 
Current formulations of the effects of a possible VLI on high energy particle 
interactions relevant in the context of UHECR, therefore, adopt a phenomeno- 
logical approach in which the form of the possible VLI is parametrized in var- 
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ious ways. VLI generally implies the existence of a universal preferred frame 
which is usually identified with the frame that is comoving with the expansion 
of the Universe, in which the CMB is isotropic. 

A direct way of introducing VLI is through a modification of the standard 
dispersion relation, E 2 — p 2 = m 2 , between energy E and momentum p = \p\ 
of particles, m being the invariant mass of the particle. Currently there is no 
unique way of parameterizing the possible modification of this relation in a 
Lorentz non-invariant theory. We discuss here a parametrization of the mod- 
ified dispersion relation which covers most of the qualitative cases discussed 
in the literature and, for certain parameter values, allows to completely evade 
the GZK limit, 

^- m ^- 2 ^-^-C^. (6) 



Here, the Planck mass Mpi characterizes non-renormalizable effects with di- 
mensionless coefficients £ and (, and the dimensionless constant d exempli- 
fies VLI effects due to renormalizable terms in the Lagrangian. The standard 
Lorentz invariant dispersion relation is recovered in the limit ~^ 0- 111 

critical string theory, effects second order in M^ 1 , ( 7^ 0, can be induced due 
to quantum gravity effects. The constants d 7^ can break Lorentz invariance 
spontaneously when certain Lorentz tensors have couplings to fermions 
of the form d^ u ifj'-f fl d u ifj, and acquire vacuum expectation values of the form 
(duv) = d8°8l [79]. Effects of first order in M^j 1 , £ 7^ 0, are possible, for 
example, in non-critical Liouville string theory due to recoiling D-branes [80]. 

Now, consider the GZK photo-pion production process in which a nucleon of 
energy E, momentum p and mass collides head-on with a CMB photon 
of energy e producing a pion and a recoiling nucleon. The threshold initial 
momentum of the nucleon for this process according to standard Lorentz in- 
variant kinematics is 

Pth,o = {ml + 2m 7r m JV )/4e , (7) 



where m n and are the pion and nucleon masses, respectively. Assuming ex- 
act energy-momentum conservation but using the modified dispersion relation 
given above, in the ultra-relativistic regime m <C p <C M, and neglecting sub- 
leading terms, the new nucleon threshold momentum p t h under the modified 
dispersion relation Eq. (6) for d — satisfies [81] 

-f3x A - ax 3 + x - 1 = , (8) 



where x = p t h/Pth,o, and 
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(ml + 2m 7T m N )M P \ (m n + m N ) 2 

2(ml + 2m 7r m Af )Mp 1 (m^ + m N ) 2 ' 

One can show that the same modified dispersion relation Eq. (6) leads to the 
same condition Eq. (8) for absorption of high energy gamma rays through 
e + e~ pair production on the infrared, microwave or radio backgrounds, if one 
substitutes p t h,o = "^/e, ot = £p?h,o/ '(8mgM P i), (3 = ?>^p^ Q / (\§m 2 e M^i where 
m e is the electron mass. 

If £, ( ~ 1, there is no real positive solution of Eq. (8), implying that the 
GZK process does not take place and consequently the GZK cutoff effect 
disappears completely. Thus UHE nucleons and/or photons will be able to 
reach Earth from any distance. On the other hand, if future UHECR data 
confirm the presence of a GZK cutoff at some energy then that would imply 
upper limits on the couplings £ and (, thus probing specific Lorentz non- 
invariant theories. If p t h — Pth,0) one could conclude from Eq. (8) that a, j3 < 1, 
which translates into |£| < 1CT 13 for the first order effects, and \(\ < 1CT 6 for 
the second order effects, £ = [81]. Confirmation of a cut-off for TeV photons 
with next-generation 7— ray observatories would lead to somewhat weaker 
constraints [82]. 

In addition, the non-renormalizable terms in the dispersion relation Eq. (6) 
imply a change in the group velocity which for the first-order term leads to 
time delays over distances r given by 

At~£r-^~n T —— J f^-) sec. (10) 



Mi 



pi 



For |£| ~ 1 such time delays could be measurable, for example, by fitting the 
arrival times of 7— rays arriving from 7— ray bursts to the predicted energy 
dependence. 

We mention that if VLI is due to modification of the space-time structure 
expected in some theories of quantum gravity, for example, then the strict 
energy-momentum conservation assumed in the above discussion, which re- 
quires space-time translation invariance, is not guaranteed in general, and 
then the calculation of the modified particle interaction thresholds becomes 
highly non-trivial and non-obvious. Also, it is possible that a Lorentz non- 
invariant theory while giving a modified dispersion relation also imposes ad- 
ditional kinematic structures such as a modified law of addition of momenta. 
Indeed, Ref. [77] gives an example of a so-called ^-Minkowski non-commutative 
space-time in which the modified dispersion relation has the same form as in 
Eq. (6) but there is also a modified momentum addition rule which compen- 
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sates for the effect of the modified dispersion relation on the particle interac- 
tion thresholds discussed above leaving the threshold momentum unaffected 
and consequently the GZK problem unsolved. In scenarios where the relativity 
of inertial frames is preserved by a non-linear representation of the Poincare 
group, thresholds are in general significantly modified only if the effective 
mass scale Mp\/xi is of the order of the unmodified threshold energy in the 
laboratory frame [78]. 

VLI by dimensionless terms such as d in Eq. (6) has been considered in 
Ref. [76]. These terms are obtained by adding renormalizable terms that 
break Lorentz invariance to the Standard Model Lagrangian. The dimension- 
less terms can be interpreted as a change of the maximal particle velocity 
Vuiax = dE/dp\E, P ?pm — 1 — d. At a fixed energy E one has the correspon- 
dence d -> (i/2)(E/M Pl ) + (C/2)(£/M P1 ) 2 , as can be seen from Eq. (6). The 
above values for £ and ( influencing the GZK effect then translate into values 
\d\ > lCT 24 . 

There are several other fascinating effects of allowing a small VLI, some of 
which are relevant for the question of origin and propagation of UHECR, 
and the resulting constraints on VLI parameters from cosmic ray observations 
are often more stringent than the corresponding laboratory limits; for more 
details, see Ref. [76] and [83]. 



6 Top-Down Scenarios 

The difficulties of bottom-up acceleration scenarios discussed earlier motivated 
the proposal of the "top-down" scenarios, where UHECRs, instead of being 
accelerated, are the decay products of certain "X" particles of mass close to the 
GUT scale. Such particles can be produced in basically two ways: If they are 
very short lived, as usually expected in many GUTs, they have to be produced 
continuously. The only way this can be achieved is by emission from TDs left 
over from cosmological phase transitions that may have occurred in the early 
Universe at temperatures close to the GUT scale, possibly during reheating 
after inflation. TDs necessarily occur between regions that are causally discon- 
nected, such that the orientation of the order parameter associated with the 
phase transition, can not be communicated between these regions and con- 
sequently will adopt different values. Examples are cosmic strings, magnetic 
monopoles, and domain walls. The defect density is thus given by the particle 
horizon in the early Universe. The defects are topologically stable, but time 
dependent motion leads to the emission of particles of mass comparable to the 
temperature at which the phase transition took place. The associated phase 
transition can also occur during reheating after inflation. 
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Alternatively, instead of being released from TDs, X particles may have been 
produced directly in the early Universe and, due to some unknown symmetries, 
have a very long lifetime comparable to the age of the Universe. In contrast 
to Weakly-Interacting Massive Particles (WIMPS) below a few hundred TeV 
which are the usual dark matter candidates motivated by, for example, su- 
persymmetry and can be produced by thermal freeze out, such super-heavy 
X particles have to be produced non-thermally (see Ref. [84] for a review). In 
all these cases, such particles, also called "WIMPZILLAs", would contribute 
to the dark matter and their decays could still contribute to UHECR fluxes 
today, with an anisotropy pattern that reflects the dark matter distribution 
in the halo of our Galaxy. However, these scenarios predict the absence of the 
GZK cutoff [85] and thus will be ruled out if the existence of the GZK cutoff 
is confirmed [8]. 

It is interesting to note that one of the prime motivations of the inflationary 
paradigm was to dilute excessive production of "dangerous relics" such as TDs 
and super-heavy stable particles. However, such objects can be produced right 
after inflation during reheating in cosmologically interesting abundances, and 
with a mass scale roughly given by the inflationary scale which in turn is fixed 
by the CMB anisotropies to ~ 10 13 GeV [84]. The reader will realize that this 
mass scale is somewhat above the highest energies observed in CRs, which 
implies that the decay products of these primordial relics could well have 
something to do with UHECRs which therefore can probe such scenarios! 

For dimensional reasons the spatially averaged X particle injection rate can 
only depend on the mass scale m x and on cosmic time t in the combination 

n x (t) = Km p x r* +P , (11) 

where k and p are dimensionless constants whose value depend on the spe- 
cific top-down scenario [86]. For example, the case p — 1 is representative of 
scenarios involving release of X particles from TDs, such as ordinary cosmic 
strings [87], monopoles connected by strings, so-called "necklaces" [88] and 
magnetic monopoles [89]. This can be easily seen as follows: The energy den- 
sity p s in a network of defects has to scale roughly as the critical density, 
p s oc p cr it oc t~ 2 , where t is cosmic time, otherwise the defects would either 
start to over-close the Universe, or end up having a negligible contribution to 
the total energy density. In order to maintain this scaling, the defect network 
has to release energy with a rate given by p s = —ap s /t oc t~ 3 , where a = 1 in 
the radiation dominated era, and a = 2/3 during matter domination. If most 
of this energy goes into emission of X particles, then typically k ~ 0(1). In the 
numerical simulations presented below, it was assumed that the X particles 
are non-relativistic at decay. 

The X particles could be gauge bosons, Higgs bosons, super-heavy fermions, 
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etc. depending on the specific GUT. They would have a mass mx comparable 
to the symmetry breaking scale and would decay into leptons and / or quarks of 
roughly comparable energy. The quarks interact strongly and hadronize into 
nucleons (Ns) and pions, the latter decaying in turn into 7-rays, electrons, and 
neutrinos. Given the X particle production rate, drix/dt, the effective injection 
spectrum of particle species a (a = 7, N, e ± , v) via the hadronic channel can 
be written as (dnx / 'dt)(2/mx)(dN a / dx), where x = 2E/mx, and dN a /dx is 
the relevant fragmentation function (FF). 

The FFs are governed by quantum chromo dynamics (QCD). At energies 
reachable at accelerators the FFs have been measured in some detail. At the 
GUT scale energies relevant for top-down scenarios, however, extrapolations 
based on QCD have to be employed to obtain the FFs. In our calculations 
we adopt the Local Parton Hadron Duality (LPHD) approximation [90] ac- 
cording to which the total hadronic FF, dN h /dx, is taken to be proportional 
to the spectrum of the partons (quarks/gluons) in the parton cascade (which 
is initiated by the quark through perturbative QCD processes) after evolving 
the parton cascade to a stage where the typical transverse momentum trans- 
fer in the QCD cascading processes has come down to ~ R^ 1 ~ few hundred 
MeV, where R is a typical hadron size. The parton spectrum is obtained from 
solutions of the standard QCD evolution equations in modified leading loga- 
rithmic approximation (MLLA) which provides good fits to accelerator data 
at LEP energies [90]. Within the LPHD hypothesis, the pions and nucleons 
after hadronization have essentially the same spectrum. The LPHD does not, 
however, fix the relative abundance of pions and nucleons after hadroniza- 
tion. Motivated by accelerator data, we assume the nucleon content /jv of 
the hadrons to be ~ 10%, and the rest pions distributed equally among the 
three charge states. Recent work on FFs goes beyond MLLA by solving the 
QCD evolution equations numerically and also includes the supersymmetric 
degrees of freedom at energies above the supersymmetry breaking scale [91]. 
The main difference in the results seems to be a nucleon-to-pion ratio that 
can be significantly higher at large x values for the extremely high energies 
of interest here [91]. However, the situation is not completely settled yet and 
the flux predictions are not very sensitive to such differences in the injection 
spectrum. 



6.1 Predicted Fluxes 

Fig. 7 shows results for the time averaged nucleon, 7— ray, and neutrino fluxes 
in a typical top-down scenario, along with constraints on diffuse 7— ray fluxes 
at GeV energies and neutrino flux sensitivities of future experiments. The 
spectrum was optimally normalized to allow for an explanation of the ob- 
served UHECR events, assuming their consistency with a nucleon or 7— ray 
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Fig. 7. [59] Predictions for the differential fluxes of 7— rays (dashed line), nucle- 
ons (dotted line), and neutrinos per flavor (solid line, assuming maximal mixing 
among all flavors) in a top-down model characterized by p = 1, rax = 2 x 10 14 GeV, 
and the decay mode X — > q + q, assuming the FF in MLLA without supersymme- 
try [90], with a fraction of 10% nucleons. The calculation used the code described in 
Ref. [28] and assumed the minimal URB version consistent with observations [36] 
and an EGMF of 10~ 12 G. Measured CR, neutrino, and 7— ray data are as in Fig. 6 
(however, only AGASA UHECR data are shown). Also shown are expected sen- 
sitivities of the currently being constructed Auger project to electron/muon and 
tau-neutrinos [92], and the planned projects telescope array (TA) [93], the fluores- 
cence/Cerenkov detector MOUNT [94], and the space based OWL [95] (we take 
the latter as representative also for EUSO [96]), the water-based NT200+ [63], 
ANTARES [97], and the ice-based ICECUBE [98], as indicated. 

primary. The flux below < 2 x 10 19 eV is presumably due to conventional 
acceleration in astrophysical sources and was not fit. The PP process on the 
CMB depletes the photon flux above 100 TeV, and the same process on the 
IR/O background causes depletion of the photon flux in the range 100 GeV- 
100 TeV, recycling the absorbed fluxes to energies below 100 GeV through 
EM cascading. The predicted background is not very sensitive to the specific 
IR/O background model, however [99]. The scenario in Fig. 7 obviously obeys 
all observational constraints within the normalization ambiguities. Note that 
the diffuse 7— ray background measured by EGRET [33] up to 10 GeV puts a 
strong constraint on these scenarios, especially if there is already a significant 
contribution to this background from conventional sources such as unresolved 
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7— ray blazars [100]. However, this constraint is much weaker for TDs or de- 
caying long lived X particles with a non-uniform clustered density [101]. 

The energy loss and absorption lengths for UHE nucleons and photons are 
< 100 Mpc. Thus, predicted UHE nucleon and photon fluxes are independent 
of cosmological evolution. The 7— ray flux below ~ 10 11 eV, however, scales as 
the total X particle energy release integrated over all redshifts and increases 
with decreasing p [102]. For m x = 2 x 10 16 GeV, scenarios with p < 1 are 
therefore ruled out, whereas models with a comovingly constant injection rate 
(p = 2) are well within the limits. 

It is clear from the above discussions that the predicted particle fluxes in the 
top-down scenarios are currently uncertain to a large extent due to parti- 
cle physics uncertainties (e.g., mass and decay modes of the X particles, the 
quark FF, the nucleon fraction f N , and so on) as well as astrophysical uncer- 
tainties (e.g., strengths of the radio and infrared backgrounds, extragalactic 
magnetic fields, etc.). More details on the dependence of the predicted UHE 
particle spectra and composition on these particle physics and astrophysical 
uncertainties are contained in Ref. [103]. 

We stress here that there are viable top-down scenarios which predict nucleon 
fluxes comparable to or even higher than the 7— ray flux at all energies, even 
though 7— rays dominate at injection. This occurs, e.g., in the case of high 
URB and/or for a strong EGMF, and a nucleon fraction of ~ 10%. Some of 
these top-down scenarios would therefore remain viable even if UHECR in- 
duced EAS should be proven inconsistent with photon primaries [104]). This 
is in contrast to scenarios with decaying massive dark matter in the Galactic 
halo which, due to the lack of interactions of injected particles, predict com- 
positions directly given by the FFs, i.e. domination by 7— rays, and thus may 
be in conflict with observed compositions [104]. 

The normalization procedure to the UHECR flux described above imposes the 
constraint Quhecr ~ 10~ 22 eV cm~ 3 sec -1 within a factor of a few [105,103,106] 
for the total energy release rate Qq from TDs at the current epoch. In most 
top-down models, because of the unknown values of the parameters involved, 
it is currently not possible to calculate the exact value of Qq from first prin- 
ciples, although it has been shown that the required values of Qo (in order to 
explain the UHECR flux) mentioned above are quite possible for certain kinds 
of TDs. Some cosmic string simulations and the necklace scenario suggest that 
defects may lose most of their energy in the form of X particles and estimates 
of this rate have been given [107,88]. If that is the case, the constraint on 
Quhecr translates via Eq. (11) into a limit on the symmetry breaking scale 
i] and hence on the mass mx of the X particle: rj ~ mx 10 13 GeV [108]. 
Independently of whether or not this scenario explains UHECR, the EGRET 
measurement of the diffuse GeV 7— ray background leads to a similar bound, 
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Qem ~ 2.2 x 10~ 23 h(3p — 1) eVcm -3 sec -1 , which leaves the bound on i] 
and mx practically unchanged. Furthermore, constraints from limits on CMB 
distortions and light element abundances from 4 He-photo-disintegration are 
comparable to the bound from the directly observed diffuse GeV 7-rays [102]. 
That these crude normalizations lead to values of rj in the right range suggests 
that defect models require less fine tuning than decay rates in scenarios of 
metastable massive dark matter. Furthermore, extragalactic top-down models 
would probably be less motivated, but not be ruled out automatically if the 
existence of the GZK cutoff is confirmed [8], in contrast to decaying massive 
dark matter scenarios. 

As discussed above, in top-down scenarios most of the energy is released in 
the form of EM particles and neutrinos. If the X particles decay into a quark 
and a lepton, the quark hadronizes mostly into pions and the ratio of energy 
release into the neutrino versus EM channel is r ~ 0.3. The energy fluence in 
neutrinos and 7— rays is thus comparable. However, whereas the photons are 
recycled down to the GeV range where their flux is constrained by the EGRET 
measurement, the neutrino flux is practically not changed during propagation 
and thus reflects the injection spectrum. Its predicted level is consistent with 
all existing upper limits (compare Fig. 7 with Fig. 6) but should be detectable 
by several experiments under construction or in the proposal stage (see Fig. 7). 
This would allow to directly see the quark fragmentation spectrum. 



7 Experimental Projects and the Future 

The field of ultra-high energy cosmic radiation is certainly experimentally 
driven. This is obvious from the lack of convergence of theoretical models 
caused by the sparseness of yet available data. To demonstrate the prospects 
for a considerable increase of data we therefore conclude with a short discussion 
of experimental projects. 

An upscaled version of the old Fly's Eye experiment, the High Resolution 
Fly's Eye (HiRes) detector is looking for CRs in Utah, USA [8]. Taking into 
account a duty cycle of about 10% because a fluorescence detector requires 
clear, moonless nights, this instrument will collect events above 10 17 eV at a 
rate about 10 times larger than the old Fly's Eye, corresponding to a few 
events above 10 20 eV per year. The largest project presently under construc- 
tion is the Pierre Auger Giant Array Observatory [21] planned for two sites, 
one in Mendoza, Argentina and another in Utah, USA for maximal sky cov- 
erage. Each site will have a ~ 3000 km ground array. The southern site will 
have about 1600 particle detectors (separated by 1.5 km each) overlooked by 
four fluorescence detectors. The ground arrays will have a duty cycle of nearly 
100%, leading to detection rates about 30 times as large as for the AGASA 
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array, i.e. about 50 events per year above 10 20 eV. About 10% of the events 
will be detected simultaneously by the ground array and the fluorescence com- 
ponent and can be used for cross calibration and detailed EAS studies. The 
detection energy threshold will be around 10 18 eV. These instruments will also 
have considerable sensitivity to neutrinos above ~ 10 19 eV, typically from the 
near-horizontal air-showers that are produced by them [92] , as shown in Fig. 7. 
The old Fly's Eye experiment [65] and the AGASA experiment [64] have al- 
ready established upper limits on neutrino fluxes based on the non-observation 
of horizontal air showers, see Fig. 6. 

As Fig. 7 demonstrates for the top-down models, neutrino fluxes in the in- 
termediate energy range from ~ 10 12 eV to ~ 10 15 eV can also are relevant 
for certain models including active galactic nuclei as CR sources. Tradition- 
ally, neutrinos in this range are detected via the EM showers or the Cerenkov 
light induced by the leptons produced by charged current reactions in ice 
or water. The largest operative experiment using ice is the Antarctic Muon 
And Neutrino Detector Array (AMANDA) [109] at the South pole which also 
has a next generation version named ICECUBE [98] which may reach up 
to ~ 10 18 eV. A water based version of this technique is operative in Lake 
Baikal [110] and another experiment for Astronomy with a Neutrino Tele- 
scope and Abyss environmental RESarch (ANTARES) is under construction 
in the Mediterranean [97]. 

In addition, there are plans to construct telescopes to detect fluorescence and 
Cerenkov light from near-horizontal showers produced in mountain targets 
by neutrinos in the intermediate window of energies between ~ 10 15 eV and 
~ 10 19 eV [111,94]. The alternative of detecting neutrino and CR induced 
showers by triggering onto the radio pulses emitted by them is also investi- 
gated currently [112]. Two implementations of this technique, the Radio Ice 
Cerenkov Experiment (RICE), a small array of radio antennas in the South 
pole ice [66], and the Goldstone Lunar Ultra-high energy neutrino Experiment 
(GLUE), based on monitoring of the moons rim with the NASA Goldstone 
radio telescope for radio pulses from neutrino-induced showers [67], have so 
far produced neutrino flux upper limits. It may be possible to use future large 
scale projects for radio telescopes such as the LOw Frequency ARray (LO- 
FAR) for UHECR and neutrino studies of low energy and angular resolution 
but high statistics [113]. Acoustic detection of neutrino induced interactions 
is also being considered [114]. 

There are also plans to detect EAS in the Earth's atmosphere from space. 
This would provide an increase by another factor ~ 50 in collecting power 
compared to the Pierre Auger Project, i.e. an event rate above 10 20 eV of up 
to a few thousand per year. Two concepts are currently being studied, the 
Orbiting Wide-angle Light-collector (OWL) [95] in the USA and the Extreme 
Universe Space Observatory (EUSO) [96] in Europe of which a prototype may 
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fly on the International Space Station. 

Space-based detectors would be especially suitable for detection of very small 
event rates such as those caused by neutrino primaries which rarely interact in 
the atmosphere due to their small interaction cross sections. This disadvantage 
for the detection process is at the same time a blessing because it makes these 
elusive particles reach us unattenuated over cosmological distances and from 
very dense environments where all other particles (except gravitational waves) 
would be absorbed. Giving rise to showers typically starting deep within the 
atmosphere, they can also be distinguished from other primaries. 

We believe that these experimental developments will allow to test, constrain, 
or rule out many of the theoretical scenarios and speculations discussed here. 
This will begin already within the next few years. 
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